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REFLECTION  OF  INCIDENT  SH  WAVE  AT  STRESS-FREE  PLANE 
BOUNDARY  TN  SFMT.TNFTNTTF  TRANSVERSELY  ISOTROPIC  MEDIUM 
WITH  PLANE  OF  ISOTROPY  PARALLEL  TO  PLANE  BOUNDARY 

1.  REFLECTED  SH  WAVE 

A  plane  progressive  SH  wave  may  be  represented  as 

(u,v,w)=A (/>„ Py,  P,)  exp{i0)(5^  +  S^y  +  S,z  - 1)}  (1) 

for  example,  see  [3],  where  m,  v  and  w  are  the  displacement  components  of  a  point  in  the 
medium  along  the  x,  y  and  z  axes,  respectively;  A  is  the  amplitude  of  particle 
displacement;  P,,  P^  and  P,  are  the  components  of  the  unit  vector  of  particle  displacement 
along  the  x,  y  and  z  axes,  respectively;  i  ■  w  denotes  radian  frequency;  S,,  and  S, 
are  the  components  of  the  slowness  vector,  which  points  in  the  same  direction  as  the 
normal  to  the  wave  front  and  has  magnitude  equal  to  the  reciprocal  of  the  magnitude  of 
the  phase  velocity,  see  [4],  along  the  x,  y  and  z  axes,  respectively;  and  t  denotes  time. 

A  plane  progressive  SH  wave  is  incident  on  the  plane  boundary  of  a  semi-infinite 
linearly  elastic  transversely  isotropic  continuum  with  plane  of  isotropy  parallel  to  the 
plane  boundary.  Define  a  cartesian  coordinate  system  (x,  y,  z)  as  follows:  the  plane 
boundary  of  the  medium  contains  the  x  and  y  axes,  and  the  z  axis  is  the  zonal  axis  of  the 
medium  (see  Fig.  1).  The  generalized  Hooke’s  Law,  when  written  relative  to  the  (x,  y,  z) 
coordinate  system  is  [4] 


+Ci2V,>’  +  Ci3W,Z 
'^yy  Cj2W^  +  CjjV  4"  CljW  ,Z 
=  Ci3«^  +  C13V  J  +  C33W,z 
'^^  =  C44(U,Z  +W^) 

V  =  ^44(V,Z+W,>’) 

'C;^  =  C66(M.3'+V^)  (2) 

where  for  i  =j  t^.  is  a  normal  stress  and  for  2  x..  is  a  shear  stress;  denotes  partial 

differentiation  with  respect  to  the  succeeding  variable;  C„,  C33,  and  are  the 

five  independent  elastic  constants  for  a  linearly  elastic  transversely  isotropic  medium; 
and  Cgg  is  equal  to  ly2(C,j-C,J. 

The  stresses  associated  with  a  plane  progressive  SH  wave  can  then  be  evaluated  by 
substituting  the  expression  for  the  displacement  components  at  a  point  in  the  medium, 
given  by  eqn.  (1),  into  eqn.  (2),  as 

+  C, 253,^3,  +  Ci3S,/’Jexp{z (0(5^  +  S^y  +  S^z  -  r)} 

X33,  —  2  cM  (Ci2‘S'j/’ ^  +  C^SyPy  +  exp{/  o)(S^  +  Syy  +  S^z  —  r)} 

—  i  coA  (C,35j/* ^  +  C^jSyPy  +  C-y^S^P j)  exp{2  co(5^  +  S^y  +  S^z  —  r)} 

=  2  ca4  {C^Sf,  +  C^SJ>,)  exp{2  (0(5^  +  +  5,z  -  r)} 

13^  =  2  0V1  {C^^Sfy  +  C^SyP^)  exp{2CO(5,^  +  S^y  +  5,z  -  r)} 

T3,  =  i(£AiCf^yP^  +  C^J^y) exp{2(o(5^  +  S^y  +  S^z  - r)}  (3) 

The  stress  boundary  conditions  on  a  stress-free  plane  boundary  require  [3] 
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<>  +  <^  =  0 

T^4<f  =  0  (4) 

where  the  represent  the  stresses  associated  with  the  incident  wave  and  the 
represent  the  stresses  associated  with  the  reflected  waves.  The  stresses  determined  by 
eqn.  (3)  satisfy  the  boundary  conditions,  eqn.  (4),  only  when  the  frequency,  co,  of  any 
reflected  wave  is  equal  to  the  frequency  of  the  incident  wave  [3]  and 

=  (5) 

where  Sj'>  and  Sf>  are  the  x  and  y  components  of  the  slowness  vector  of  the  incident 
wave,  respectively,  and  Sf>  and  are  the  x  and  y  components  of  the  slowness  vector  of 
any  reflected  wave,  respectively  [2],  Eqn.  (5)  indicates  that  the  incident  and  the  reflected 
waves  lie  in  the  same  plane,  called  the  plane  of  incidence.  Without  loss  of  generality,  the 
plane  of  incidence  can  be  assumed  to  coincide  with  the  y-z  plane;  the  analysis  is  thus 
simplified.  Then,  the  x  components  of  the  slowness  vectors  of  the  incident  and  reflected 
waves  can  be  allowed  to  vanish,  thereby  satisfying  the  relation 

5«  =  5f>  =  0  (6) 

An  SH  wave  with  slowness  surface  in  a  plane  containing  the  zonal  axis  of  a 
transversely  isotropic  medium  possesses  a  transverse  displacement  [5].  Thus,  for  an  SH 
wave  traveling  in  the  y-z  plane  in  the  (x,  y,  z)  coordinate  system,  the  unit  vector  of 
particle  displacement  is  given  by 


(/>„/’,,/>J  =  (l,0,0) 


5 


A  wave  with  P^,  and  P^  given  by  eqn.  (7)  and  x  component  of  slowness  vector  given  by 
eqn.  (6)  has  stresses  determined  from  eqn.  (3)  satisfying 

X  =x  =1  =T  =0 

''xx  ^yy  ^zz  ^yz  ^ 

(8) 

A  P  wave  with  slowness  surface  in  a  plane  containing  the  zonal  axis  of  a  transversely 
isotropic  medium  possesses  a  vanishing  transverse  displacement  [5].  Therefore,  a  P  wave 
traveling  in  the  plane  x=0  in  the  (x,  y,  z)  coordinate  system  has  a  unit  vector  of  particle 
displacement 


(P^,P^,P,)  =  (0,Pp,Pf)  (9) 

where  Pf>  and  Pj'’>  are,  respectively,  the  directionally  dependent  y  and  z  components  of 
the  unit  vector  of  particle  displacement  of  a  P  wave. 

A  wave  P^,  P^  and  P  satisfying  eqn.  (9)  and  determined  by  eqn.  (6)  has  stresses 
given  by  eqn.  (3)  satisfying 


(10) 

An  SV  wave  with  slowness  surface  in  a  plane  containing  the  zonal  axis  of  a 
transversely  isotropic  medium  also  has  a  vanishing  transverse  displacement  [5].  The  unit 
vector  of  particle  displacement  of  an  SV  wave  traveling  in  the  y-z  plane  in  the  (x,  y,  z) 
coordinate  system  is  then 


(11) 
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where  and  are,  respectively,  the  directionally  dependent  y  and  z  components  of 
the  unit  vector  of  particle  displacement  of  an  S  V  wave. 

A  wave  with  unit  vector  of  particle  displacement  determined  by  eqn.  (11)  and  x 
component  of  slowness  vector  given  eqn.  (6)  has  stresses  evaluated  from  eqn.  (3) 
satisfying 


(12) 


An  SH  wave  is  incident  on  the  plane  boundary  of  a  semi-infinite  transversely  isotropic 
medium.  The  stresses  at  the  point  of  incidence  associated  with  the  incident  SH  wave  are 
constrained  by  eqn.  (8)  and  satisfy 


^yz  ^zz  ^ 

(13) 

where  the  superscript  ^  is  used  to  denote  properties  associated  with  the  incident  wave. 


The  stress  boundary  conditions  for  any  wave  reflected  by  an  incident  SH  wave  can  be 
expressed  by  using  eqn.  (4)  and  eqn.  (8)  as 


=  (14) 

where  the  superscript  «>  is  used  to  denote  properties  associated  with  the  reflected  wave. 
Eqn.  (10)  indicates  that  the  Xy,  and  components  of  the  stress  tensor  associated  with  a  P 
wave  are  not  zero  and  eqn.  (12)  indicates  that  the  Xy,  and  x„  components  of  the  stress 
tensor  associated  with  an  SV  wave  are  also  nonzero.  Therefore,  no  P  wave  or  SV  wave 
will  be  reflected  back  into  the  medium  by  an  incident  SH  wave,  because  a  reflected  P 
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wave  or  SV  wave  would  have  associated  stresses  that  would  violate  the  stress  boundary 
conditions.  Eqn.  (8)  indicates  that  a  reflected  SH  wave  can  satisfy  the  boundary 
conditions.  Therefore,  an  incident  SH  wave  will  cause  only  an  SH  wave  to  be  reflected 
back  into  the  medium,  as  shown  in  [2]. 


2.  SLOWNESS  SURFACE  OF  SH  WAVE 

The  equations  of  motion  relative  to  the  (x,  y,  z)  coordinate  system  are  [4] 

+'^yz>y  +  =  pw,tt  (15) 

where  p  is  the  density  of  the  medium  and  ",tt"  denotes  the  second  derivative  of  the 
preceeding  variable  with  respect  to  time. 

When  the  components  of  stress  are  calculated  according  to  eqn.  (3)  and  the 
displacement  components  are  calculated  according  to  eqn.  (1),  the  equations  of  motion 
can  be  written  as  [2] 


(CnS^  +  C^y  +  -  9)P,  +  {C,,  +  C  +  {C,,  +  =  0 

(Ci2  +  C^)S.S^P,  +  (C^S^  +  -  p)P^  +  (C,3  +  =  0 

(Ci3  +  +  (C,3  +  CJS^S^P^  +  {C^sl  +  c^sl  +  C33S,"  -  p)F,  =  0 

Eqn.  (16)  can  be  rewritten  in  matrix  form  as 


where  [5]  is  a  3x3  matrix  with  entries 


(17) 
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bn  =  iCnSl  +  +  C^Sl  -  p) 

b^^  =  {C^l  +  CnSl+C^S^.-9) 

^33  =  iC^^l  +  <^44^^  +  C^^Sl  -  p) 
bn  —  ^21  “  (^12  ■*■  ^66)^x^y 
bn  ~  bn  =  (Cj3  +  C44)S,5j 

&23  =  ^32  =  (Cl3  +  CJV,  (18) 

{P  JP)  is  the  unit  vector  of  particle  displacement  and  [0]  is  the  3x1  zero  matrix.  The 
matrix  [fi]  as  defined  by  eqn.  (18)  is  a  symmetric  matrix;  therefore,  there  exist  three  real 
eigenvalues.  One  eigenvalue  corresponds  to  an  SH  wave,  another  eigenvalue  corresponds 
to  a  P  wave,  and  the  third  eigenvalue  corresponds  to  an  SV  wave  [2]. 

The  plane  wave  solution  is  found  by  setting  the  determinant  of  matrix  [fi],  the  matrix 
of  coefficients  of  P,,  P,  and  P,,  equal  to  zero  [6].  Expanding  the  determinant  of  [B]  and 
solving  for  the  three  roots  yields  equations  for  the  three  slowness  surfaces.  The  slowness 
surface  for  an  SH  wave  is  [7] 

C^(Sl  +  S^)  +  C^S^  =  p  (19) 

The  values  of  the  material  constants  C«,  and  p  are  given  for  a  representative 
fiberglass  epoxy  composite  in  [4]  as  C(j=3. 243x10®  N/m*,  C^=4.422xl0®  N/m^  and  p 
=1850  kg/m®.  The  slowness  surface  of  an  SH  wave  traveling  in  the  y-z  plane  in  the 
fiberglass  epoxy  composite  is  obtained  by  substituting  the  numerical  values  of  the 
constants  into  eqn.  (19)  and  setting  the  x  component  of  the  slowness  vector  equal  to  zero. 
The  first  quadrant  of  the  slowness  surface  thus  obtained  is  shown  in  Fig.  2. 
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3.  ANGLE  OF  REFLECTION  OF  REFLECTED  SH  WAVE 

According  to  eqn.  (5)  the  y  component,  S^<'> ,  of  the  slowness  vector  of  the  incident  SH 
wave  is  equal  to  the  y  component,  ,  of  the  slowness  vector  of  the  reflected  SH  wave. 
The  z  component,  5  w  ,  of  the  slowness  vector  of  the  reflected  SH  wave  is  found  from  the 
equation  for  the  slowness  surface  of  an  SH  wave,  eqn.  (19),  by  setting  equal  to 
and  Sf>  equal  to  zero,  according  to  eqn.  (6).  Then,  S<^>  is  found  to  satisfy 

(20) 

where  Sf>  is  the  z  component  of  the  slowness  vector  of  the  incident  SH  wave.  The 
slowness  vector  of  the  incident  wave  points  out  of  the  medium,  and  the  slowness  vector 
of  the  reflected  wave  points  into  the  medium;  therefore,  the  relationship  between  S  w  and 
Sj'>  must  be  [2] 


(21) 


The  angle  of  incidence  is  defined  as  the  angle  between  the  slowness  vector  of  the 


incident  wave  and  the  normal  to  the  boundary  at  the  point  of  incidence.  The  angle  of 
incidence  0^'^  of  the  incident  SH  wave  is  then 


9^^^  =  arctan 


(22) 


where  the  point  of  incidence  is  taken  to  be  on  the  plane  boundary  of  the  semi-infinite 
body  (see  Fig.  1). 

The  angle  of  reflection  is  defined,  in  a  manner  similar  to  the  angle  of  incidence,  as  the 
angle  between  the  slowness  vector  of  the  reflected  wave  and  the  normal  to  the  boundary 
at  the  point  of  reflection.  The  angle  of  reflection  0^*^  of  the  reflected  SH  wave  is  then 
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(23) 


=  arctan 


Evaluating  6^  using  eqn.  (5)  and  eqn.  (21)  and  comparing  the  resulting  expression  to 
eqn.  (22)  establishes 

q(R)^q(I)  (24) 

Thus,  the  angle  of  reflection  of  the  reflected  SH  wave  is  equal  to  the  angle  of  incidence  of 
the  incident  SH  wave  [2]. 


INPUT-OUTPUT  CHARACTERIZATTON  OF  FIBER  COMPOSITE 


1.  PHASE  VELOCITY  AND  DELAY  TIME 

Certain  fiber  composites  may  be  modeled  as  a  linearly  elastic  transversely  isotropic 
homogeneous  continua  [1].  A  fiber  composite  modeled  as  such  a  solid  in  the  form  of  an 
infinite  plate  where  the  plane  of  isotropy  lies  in  the  midplane  of  the  plate  is  to  be  studied. 
A  cartesian  coordinate  system  (x,  y,  z)  is  chosen  such  that  the  x-y  plane  coincides  with 
the  plane  of  isotropy  and  the  plate  is  bounded  above  by  the  plane  z  =  hH  and  bounded 
below  by  the  plane  z  =  -  /i/2. 

A  transmitting  transducer  and  a  receiving  transducer  are  assumed  to  be  coupled  to  the 
top  face  of  the  plate  in  the  y-z  plane  and  separated  by  a  distance  L.  The  input  electrical 
voltage  of  the  transmitting  transducer  is  a  known  function  of  time,  V.(r).  The  output 
electrical  voltage  of  the  receiving  transducer  is  an  unknown  function  of  time,  The 
transmitting  transducer  converts  the  input  electrical  voltage  into  a  stress  that  travels 
through  the  plate  as  stress  waves.  The  receiving  transducer  converts  the  stress  associated 
with  the  stress  waves  into  an  output  voltage.  In  the  following  analysis,  only  those  stress 
disturbances  at  the  receiving  transducer  associated  with  SH  waves  generated  by  the 
transmitting  transducer  will  be  considered. 

The  SH  waves  produced  by  the  transmitting  transducer  experience  multiple  reflections 
at  each  face  of  the  plate  before  reaching  the  receiving  transducer.  Since  the  isotropic 
plane  of  the  plate  lies  in  the  midplane  of  the  plate  and  is  parallel  to  the  top  and  bottom 
faces  of  the  plate,  each  reflection  may  be  treated  as  the  reflection  of  a  plane  progressive 
SH  wave  on  the  plane  boundary  of  a  semi-infinite  linearly  elastic  tranversely  isotropic 
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continuum  with  its  isotropic  plane  parallel  to  the  boundary  [2].  Therefore,  at  each 
reflection  the  incident  SH  wave  produces  only  a  reflected  SH  wave  and  the  angle  of 
reflection  of  the  reflected  SH  wave  is  equal  to  the  angle  of  incidence  of  the  incident  SH 
wave. 

The  transmitting  transducer  and  the  receiving  transducer  are  coupled  to  the  same  face 
of  the  plate,  therefore,  any  SH  wave  traveling  from  the  transmitting  transducer  to  the 
receiving  transducer  must  experience  an  odd  number  of  reflections.  Let  the  total  number 
of  reflections  experienced  by  an  arbitrary  SH  wave  traveling  from  the  transmitting 
transducer  to  the  receiving  transducer  be  equal  to  2n-l,  where  n  is  called  the  reflection 
index  and  is  a  positive  integer  equal  to  the  number  of  reflections  from  the  bottom  face  of 
the  plate.  Assume  the  angle  of  incidence  of  the  SH  wave  at  the  receiving  transducer  is  0,. 
The  angle  of  reflection  is  equal  to  the  angle  of  incidence  at  each  reflection  experienced 
by  the  SH  wave,  so  the  angle  of  incidence  at  each  reflection  is  0,  (see  Fig.  3).  The 
distance  traveled  by  the  SH  wave  between  each  reflection  is  calculated  from  Fig.  3  as 


(25) 


The  distance  traveled  in  the  y  direction  by  the  SH  wave  between  each  reflection  is 

i,„  =  r,sin0,  (26) 

Using  eqn.  (25)  to  evaluate  r  in  eqn.  (26),  the  distance  h,  can  be  expressed  as 

h,=Atan0,  (27) 
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The  SH  wave  travels  a  distance  in  the  y  direction  between  each  reflection,  between 
the  transmitting  transducer  and  the  initial  reflection,  and  between  the  final  reflection  and 
the  receiving  transducer,  so  the  total  distance  5,  traveled  in  the  y  direction  by  the  SH 
wave  is 


fi„=2«Atane„  (28) 

The  transmitting  transducer  and  the  receiving  transducer  are  separated  by  a  distance  L  in 
the  y  direction,  therefore;  the  total  distance  traveled  by  the  SH  wave  must  be  L.  Eqn.  (28) 
can  then  be  expressed  as 


L  =  2nh  tan  6, 

The  angle  of  incidence  at  the  receiving  transducer  must  therefore  satisfy 
0„  =  arctan 


(29) 


(30) 


where  n  is  a  positive  integer. 

Eqn.  (30)  allows  all  the  trigonometric  functions  of  0^  to  be  calculated  (see  Fig.  4).  The 
tangent  of  0^  is  given  by 


(31) 


The  sine  of  0,  is  given  by 
sin0,= 

and  the  cosine  of  0^  is  given  by 


jL^  +  4n^h^ 


(32) 
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2nh 

cos  9.  =  ,  ~ 

Vl  +4/1  /i 


The  distance  r,  traveled  by  the  SH  wave  between  each  reflection  can  be  evaluated 
using  eqn.  (25)  and  eqn.  (33)  as 


/Z^+4n^ 


The  SH  wave  travels  a  distance  between  each  reflection,  between  the  transmitting 
transducer  and  the  first  reflection,  and  between  the  last  reflection  and  the  receiving 
transducer^  so  the  total  distance  traveled  by  the  SH  wave  between  the  transmitting 
transducer  and  the  receiving  transducer  is 

+  (35) 

The  phase  velocity  C^(Q)  of  an  SH  wave  traveling  in  a  transversely  isotropic 
continuum  is  [4] 

Q„(0)  =  {C«sin"0  +  C^cos"0}‘''  (36) 

Using  eqn.  (32)  and  eqn.  (33),  the  phase  velocity  can  be  expressed  in  terms  of  n  as 


I  r  1  r  r  4/t^/t^ 

Q//(0«)=iC66  ^2  +  4„2;j2 


The  time  delay  r,  is  defined  as  the  time  taken  for  the  SH  wave  to  reach  the  receiving 
transducer  after  being  produced  by  the  transmitting  transducer.  The  time  delay  is  thus 


'■  CM 
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Eqn.  (38)  can  be  evaluated  using  eqn.  (35)  and  eqn.  (37)  as 

=  {L^  +  An\^}  +  (39) 

The  time  delay  for  a  typical  transducer  geometry  is  shown  in  Table  1.  The  transducers 
are  assumed  to  be  coupled  to  a  fiberglass  epoxy  composite  plate  having  continuum  model 
with  material  constants  C^=4.422xlO  N/m^  and  C^=3.234xlO  N/m^  and  separated  by  a 
distance  £/=10  cm.  The  plate  is  assumed  to  be  of  thickness  /?=5  cm.  The  time  delay  is 
shown  as  a  function  of  the  reflection  index  for  values  of  n  from  1  to  9. 

2.  ASSUMPTIONS  ON  TRANSDUCERS 

The  transmitting  transducer  and  the  receiving  transducer  are  assumed  to  be  transverse 
transducers  that  transform  an  electrical  voltage  into  a  uniform  shear  stress  or  a  uniform 
shear  stress  into  an  electrical  voltage,  respectively.  The  following  approach  parallels  that 
of  [8]. 

If  an  input  voltage  U  of  amplitude  V  and  frequency  co  is  applied  according  to 

U,  =  V  exp{-/tor}  (40) 

where  and  t  denotes  time,  the  shear  stress  that  is  introduced  into  the  medium  at 
the  transducer-medium  interface  by  the  transmitting  transducer  is 

'':«(0  =  E,(to)Uexp{-i(cor-i-(|),)}  (41) 
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where  Fj(a))  is  the  frequency  dependent  transduction  ratio  for  the  transmitting  transducer 
in  transforming  a  voltage  into  a  stress  and  (|)i  is  a  phase  angle.  In  eqn.  (43)  and  eqn.  (44) 
the  harmonic  character  of  the  signals  is  expressed  in  complex  notation,  but  only  the  real 
parts  of  these  and  subsequent  equations  should  be  considered.  The  amplitude  T  of  the 
applied  stress  is  then 

T  =  F,(G))V  (42) 

Similarly,  if  a  stress  wave  producing  a  shear  stress  component  of  amplitude  T '  and 
frequency  co  that  impinges  on  the  receiving  transducer  is  defined  as 

=  (43) 

then  the  output  voltage  from  the  receiving  transducer  is 

V„(0  =  F2(o))7’'exp{-i  (cor + (44) 
where  Fj((o)  is  the  frequency  dependent  transduction  ratio  for  the  receiving  transducer  in 
transforming  a  shear  stress  to  a  voltage,  and  is  a  phase  angle.  Thus,  the  amplitude  V'  of 
the  output  electrical  voltage  is 

V'  =  F2i(0)r  (45) 

The  characteristics  of  F^(d)  and  F^(oi)  are  unknown  except  that  the  prcxiuct  F^((£i)F^((o)  is 
dimensionless. 

3.  DIRECTIVITY  FUNCTION 

The  directivity  function  of  the  shear  stress,  t„,  due  to  an  SH  wave  traveling  in  a 
transversely  isotropic  medium  is  evaluated  from  the  far-field  asymptotic  solutions  of  the 
displacement  components  produced  by  a  harmonic  point  load  imbedded  in  an  infinite 
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body  in  [2].  The  directivity  function  is  given  by  [2] 


where  5/  and  5/  are  the  y  and  z  components,  respectively,  of  a  point  on  the  slowness 
surface  of  an  SH  wave  where  the  normal  to  the  slowness  surface  is  parallel  to  the  line 
connecting  the  transmitting  transducer  and  the  point  at  which  the  directivity  function  is  to 
be  evaluated. 


The  shear  stress,  radiated  by  the  transmitting  transducer,  that  is  incident  at  the 
receiving  transducer  with  angle  of  incidence  0,,  is  assumed  to  be  equivalent  to  the  shear 
stress  at  a  point  associated  with  an  SH  wave  traveling  in  a  semi-infinite  transversely 
isotropic  medium  (see  Fig.  5).  Except  for  the  reflection  coefficients  at  each  face  of  the 
plate  the  shear  stress  can  be  computed  as  if  there  were  no  bottom  boundary  to  the  plate 
[2]. 


The  slope,  m,  of  the  line  connecting  the  transmitting  transducer  to  the  point  is 
found  from  Fig.  5  to  be 


Ay  2nh 
^  ~Tz~~~L 


(47) 


Eqn.  (47)  establishes  that  the  line  connecting  the  transmitting  transducer  to  the  point 
lies  in  the  y-z  plane.  Therefore,  a  point  on  the  slowness  surface,  (5^,’,  S^/,  5,/ ),  where  the 
normal  is  parallel  to  the  line  connecting  the  transmitting  transducer  to  the  point  may 
be  found  by  setting  5,/  equal  to  zero. 
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The  normal  to  the  slowness  surface  of  an  SH  wave  is  found  at  any  point  in  the  y-z 
plane  by  setting  the  x  component  of  the  slowness  vector,  5,,  equal  to  zero  and 
differentiating  the  equation  for  the  slowness  surface,  eqn.  (19),  with  respect  to  S,,  the  z 
component  of  the  slowness  vector.  Solving  the  subsequent  equation  for  the  slope, 
-dSJdS^,  of  the  normal  yields 

dSy  C44SJ 

dSj  Cf^y 

Combining  eqn.  (47)  and  eqn.  (48)  establishes  the  constraint  on  5^/  and  S,,*,  the  point  on 
the  slowness  surface  where  the  normal  is  parallel  to  the  line  connecting  the  transmitting 
transducer  to  the  point  M,,  as 


The  point  (5^,,*,  S,,.*,  5,,*)  must  lie  on  the  slowness  surface  of  an  SH  wave  and  have 
slope  satisfying  eqn.  (49).  Therefore,  5,/  and  S,/  can  be  found  by  setting  5,/  equal  to 
zero  and  then  solving  eqn.  (19)  and  eqn.  (49)  simultaneously.  Then  the  value  of  5,,,  is 


and  the  value  of  S,*  is 


(51) 
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Substituting  eqn.  (50)  and  eqn.  (51)  into  eqn.  (46)  produces  an  expression  for 
the  value  of  the  directivity  function  associated  with  an  angle  of  incidence  0,  at  the 
receiving  transducer.  The  expression  for  as  a  function  of  n  is 


~  L  2tz 


.3/2 


Caa  +  C( 


66 


4n^h‘ 


'44\ 


1  + 


-|l/2 


(52) 


The  directivity  function  for  an  SH  wave  traveling  in  a  transversely  isotropic  plate  is 
calculated  as  a  function  of  the  angle  of  incidence  at  the  receiving  transducer  by  treating 
eqn.  (30)  and  eqn.  (52)  as  a  pair  of  parametric  equations  with  parameter  n  equal  to  the 
reflection  index.  The  directivity  function  is  shown  for  an  SH  wave  traveling  in  a 
fiberglass  epoxy  composite  plate  as  a  function  of  the  angle  of  incidence  at  the  receiving 
transducer  in  Figs.  6,  7  and  8.  The  transducer  geometry  is  assumed  to  be  L= 10  cm  and 
h=5  cm.  The  composite  plate  has  with  material  properties  C^=4.422xl0®  N/m^  and 
C«=3.234x105  N/m^  [4].  The  directivity  function  is  shown  for  frequencies  equal  to  0.75, 
1.00  and  1.25  MHz. 


4.  STRESS  FIELD  RADIATED  BY  TRANSMITTING  TRANSDUCER 

If  there  were  no  bottom  boundary  to  the  plate  and  the  SH  waves  were  propagating  in 
an  infinite  half-space,  the  SH  wave  incident  at  the  receiving  transducer  with  angle  of 
incidence  0,  would  travel  to  a  point  M,,  in  a  time  r,,  traveling  a  distance  /?„  (see  Fig.  6). 
The  amplitude  of  the  hypothetical  shear  stress  at  the  point  M,  is  and  is  defined  as  [2] 

£)(5«) 

=  (53) 
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where  T  is  the  magnitude  of  the  point  load  applied  by  the  transmitting  transducer  in  the  z 
direction,  Df”>,  the  value  of  the  directivity  function  for  an  SH  wave  associated  with  the 
point  M  is  given  by  eqn.  (52)  and  cc  is  the  SH  wave  attenuation  constant  of  the  medium. 
The  SH  wave,  however,  does  not  propagate  in  an  infinite  half-space,  but  instead 
experiences  2n-l  reflections.  The  amplitude  T,,  of  the  shear  stress  at  the  receiving 
transducer  is  thus  obtained  by  modifying  eqn.  (53)  as 


T„  = 


Rn 


exp{-a/?„} 


(54) 


where  is  the  reflection  coefficient  for  SH  waves  and  is  equal  to  -1  [2]. 

The  amplitude  of  the  output  voltage  from  the  receiving  transducer  is  obtained  from 
eqn.  (45)  and  eqn.  (54)  and  is 

-F2(co)rDf»> 


V.  =■ 


R. 


-exp{-cd?J 


(55) 


Substituting  eqn.  (42)  into  eqn.  (55)  yields 

-Fi(0))F2(0))VDr> 


=- 


R. 


■exp{-cd?,} 


(56) 


where  V  is  the  amplitude  of  the  input  voltage.  Introducing  a  possible  electrical  signal 
amplification  factor,  K,  eqn.  (56)  can  be  written  as 


V„ 


-ji:Fi(co)F,(a))VZ)r> 

Rn 


exp{-aFJ 


(57) 
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5.  OUTPUT  VOLTAGE  DUE  TO  TONE  BURST 


Assume  the  input  voltage  V,  to  the  transmitting  transducer  is  a  periodic  function  of 
time,  of  center  frequency  co  and  duration  r.  Mathematically  the  input  voltage  can  be 
expressed  as  the  sum  of  two  periodic  functions  of  frequency  co  that  are  180°  out  of  phase, 
one  beginning  at  time  r=0  and  one  beginning  at  time  r=r,.  The  input  voltage  is  then  given 
by 


V,  =  V  exp{/co/}  U(r)  -  V  exp{ia)r}  U(r  - 1^)  (58) 

where  V  is  the  amplitude  of  the  input  voltage,  ,  and  U(x)  is  the  unit  step  function 

defined  as 


fo 


U(x)  = 


if  x<0 


(59) 


U  if  x>0 

The  periodic  nature  of  the  signal  in  eqn.  (58)  and  subsequent  equations  is  expressed  in 
complex  notation,  but  only  the  real  part  of  the  signal  is  to  be  considered. 


The  voltage  output  V  at  the  receiving  transducer  is  found  by  considering  the  two 
terms  on  the  right-hand  side  of  eqn.  (61)  independently  and  superposing  the  voltage 
output  associated  with  each  input  signal.  Consider  the  input  voltage  V'  given  by  the  first 
term  on  the  right-hand  side  of  eqn.  (58) 


V/  =  Vexp{2(Dr}U(f)  (60) 

The  stress  waves  produced  by  a  transmitting  transducer  with  input  voltage  V'  will  have  a 
periodic  character  with  respect  to  time.  Therefore,  the  amplitude  of  the  output  voltage  V,; 
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associated  with  the  SH  wave  experiencing  2n-\  reflections  will  also  be  periodic  in  time 
and  the  output  voltage  V'  at  the  receiving  transducer  must  be  found  by  the  superposition 
of  the  contributions  of  each  reflection  path. 

An  SH  wave  experiencing  2n-\  reflections  between  the  transmitting  transducer  and 
the  receiving  transducer  will  be  in  phase  at  the  receiving  transducer  with  an  SH  wave 
traveling  in  a  semi-infinite  body  to  a  point  (see  Fig.  5).  The  stress  associated  with  a 
plane  progressive  SH  wave  is  given  by  eqn.  (3)  and  can  be  calculated  at  point  as 


where  is  the  amplitude  of  the  stress  and  Sy^  and  5,^  are  the  y  and  z  components  of  the 
slowness  vector  of  the  SH  wave,  respectively.  The  shear  stress  at  the  receiving 
transducer  is  then 

exp|z[© [S,^L  - S,2nh  - 1 

where  is  the  amplitude  of  the  stress  and  is  given  by  eqn.  (54).  The  output  voltage 
at  the  receiving  transducer  is  found  by  substituting  eqn.  (62)  into  eqn.  (44)  to  be 

n.'«)  =F,«f))r„exp|i[a)(^S,  L -S,2nh-ty</_ 

where  <])  is  a  phase  angle.  The  output  voltage  is  then  calculated  using  eqn.  (42)  and  eqn. 
(54)  and  including  a  possible  electrical  signal  amplification  factor  K  as 
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(64) 


,_-/rF,(o))F,(co)rar  .  _  fr  ^  ^ 

,  -  ^ - exp{-aFJexpji  (Si\^SyL  - S^lnh  -t\  +  (^ 


where  V  is  the  amplitude  of  the  input  voltage  and  the  is  given  by  eqn.  (52).  The  total 
output  voltage  V/  is  found  by  summing  eqn.  (64)  over  all  n  and  retarding  each 
contribution  by  the  delay  time  r  to  be 


K  - - - exp{-aFJx 


exp|i[co|^5,/  - S,^2nh  -  + 4>J|  U(r  - rj  (65) 

where  r,  is  given  by  eqn.  (38). 

The  output  voltage  Y"  associated  with  the  input  voltage  characterized  by  the  second 
term  on  the  right-hand  side  of  eqn.  (58)  is  found  similarly  to  be 

Y„  -  1 - - - exp{-ot/?„}x 

""=1  A. 


expj  i[o)  ^5,  L  -  S,2nh  -  r j  -h  0 j  I  u(r  -  (r,  -  rj) 


where  is  the  duration  of  the  input  signal. 

In  order  to  evaluate  the  infinite  sum  in  eqn.  (65),  5^,  and  must  be  found  as  functions 
of  n.  The  components  of  the  slowness  vector  (0,  SJ  associated  with  the  SH  wave 
experiencing  2n-l  reflections  between  the  transmitting  transducer  and  the  receiving 
transducer  are  constrained  by  eqn.  (22)  to  satisfy 
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(67) 


0„  =  arctan 

Substituting  eqn.  (30)  into  eqn.  (67)  yields 

=  ^  (68) 

Inh 

Solving  eqn.  (1 1)  and  eqn.  (68)  simultaneously  determines  5^,  and  5„  as 


The  output  voltage  V,  associated  with  the  input  voltage  determined  by  eqn.  (58)  is 
found  by  adding  eqn.  (65)  and  eqn.  (66).  Substituting  in  the  values  of  the  constants  given 
by  eqns.  (35),  (39)  and  (69)  allows  V.  to  be  written  as 
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V^{t)  =  -KF,{(xi)F^mV  I<^-p=^=exp{-aVz7+4rf}x 


expi 


CO 


Inh 


«  =  H  VL^4-4n^/2^ 

r  ^1/2  ^ 

+  2a2/z 


/2 


^66 


Ca 


\l/2 


4n2A^ 


V 


-t 


+0 


U|  t-{L^  +  An'h^}  {CJ.^  +  C^An\^} 


2,  2-, -1/2 


U^r  -  [{L^  +  4/z"/i'}  +  rj  j 


(70) 


where  D/"i  is  a  function  of  n  according  to  eqn.  (52). 

The  output  voltage  predicted  by  eqn.  (70)  is  shown  as  a  function  of  time  in  Figs.  9,  10 
and  11.  The  transducers  are  assumed  to  be  coupled  to  a  representative  fiberglass  epoxy 
composite  plate.  The  composite  plate  can  be  modeled  as  a  transversely  isotropic 
continuum  plate  with  material  properties  p=1850  kg/m^,  C^=4.422xlO  N/m^  and 
Cgg=3.234xl0«  N/m^  and  thickness  h=5  cm  and  transducer  separation  L=10  cm  [4].  The 
continuum  plate  is  assumed  to  have  SH  wave  attenuation  constant  a  equal  to  zero.  The 
normalized  output  voltage  VJKF^F^  is  shown  for  center  frequencies  equal  to  0.75,  1.00 
and  1.25  MHz.  and  a  signal  duration  of  60ps;  the  phase  angle  <})  is  set  equal  to  zero. 


26 


CONCLUSION  AND  DISCUSSION 


The  input-output  characterization  of  certain  fiberglass  epoxy  composite  plates  can  be 
studied  using  a  transversely  isotropic  continuum  plate  model  when  the  wavelengths  under 
consideration  are  long  compared  to  the  mean  fiber  diameter.  For  wavelengths  close  to  or 
less  than  the  mean  fiber  diameter,  a  continuum  plate  model  is  no  longer  appropriate  and 
the  inhomogeneities  of  the  composite  must  be  considered. 

In  the  analysis  of  a  transversely  isotropic  continuum  plate,  the  angle  of  reflection  of  a 
reflected  SH  wave  is  equal  to  the  angle  of  incidence  of  the  incident  SH  wave  and  only  an 
SH  wave  is  reflected  by  the  incident  SH  wave  for  reflections  at  the  top  or  bottom  face  of 
the  plate.  This  is  due  to  the  fact  that  the  isotropic  plane  of  the  plate  is  parallel  to  the  two 
faces  of  the  plate.  When  this  geometry  is  present,  the  plate  can  be  modeled  as  a 
semi-infinite  transversely  isotropic  medium,  except  for  the  reflection  coefficients 
encountered  at  each  reflection  which  must  be  factored  in  separately.  Treating  the 
continuum  plate  as  a  semi-infinite  body  reduces  the  complexity  of  tracing  stress  waves 
through  the  medium.  If  this  special  geometry  does  not  exist,  the  use  of  this  simplified 
analysis  is  inappropriate  and  the  plate  must  be  analyzed  in  a  more  complicated  manner. 

It  is  possible  to  compute  the  relationship  between  the  input  voltage  to  the  transmitting 
transducer  and  the  output  voltage  from  the  receiving  transducer  in  terms  of  the  material 
constants,  the  geometry  of  the  transducer  arrangement  and  a  single  variable,  the  reflection 
index  n.  Using  the  fact  that  the  angle  of  incidence  of  an  incident  SH  wave  is  equal  to  the 
angle  of  reflection  of  the  reflected  SH  wave,  the  angle  0  with  which  the  final  wave  of  any 
series  of  waves  from  the  transmitting  transducer  to  the  receiving  transducer  is  incident  at 
the  receiving  transducer  can  be  computed  as  a  function  of  the  reflection  index.  The 
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nonzero  components  of  the  slowness  vector  of  an  SH  wave  traveling  in  a  plane 
perpendicular  to  the  isotropic  plane  of  the  plate,  the  phase  velocity  of  the  SH  wave,  the 
time  delay  in  reaching  the  receiving  transducer,  the  directivity  function  of  the  shear  stress 
produced  by  the  transmitting  transducer  and  the  output  voltage  from  the  receiving 
transducer  can  then  be  computed  as  functions  of  the  reflection  index  and  the  geometric 
and  physical  properties  of  the  plate. 

The  output  voltage  from  the  receiving  transducer,  as  predicted  by  eqn.  (70),  is  shown 
in  Figs.  9,  10  and  11  as  functions  of  time  for  input  voltages  of  60  |is  duration  and  center 
frequencies  of  0.75,  1.00  and  1.25  MHz.,  respectively.  All  three  output  voltages  are 
observed  to  exhibit  similar  behavior,  being  characterized  by  simple  harmonic  functions 
experiencing  discontinuities  in  amplitude.  Checking  the  times  at  which  these  changes  in 
amplitude  occur  against  Table  1,  which  shows  the  time  delays  for  reflection  indices  from 
1  to  9,  reveals  that  the  amplitude  changes  occur  at  the  time  delays  for  these  values  of  the 
reflection  index.  Thus  the  discontinuities  in  the  output  voltages  correspond  with  the 
arrival  times  of  different  wave  paths.  Since  the  time  delay  is  independent  of  the 
frequency  of  the  wave,  these  discontinuities  occur  at  the  same  time  in  each  output  voltage 
regardless  of  the  frequency  of  the  associated  input  voltage. 

This  study  enhances  the  theoretical  understanding  of  the  nondestructive  evaluation 
(NDE)  of  transversely  isotropic  media  such  as  certain  fiber  composites. 
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Time  delay  in  reaching  receiving  transducer  for 
wave  in  composite  plate  for  reflection  indices 
from  1  to  9 . 


reflection  index 


1 

2 

3 

4 

5 

6 

7 

8 
9 


time  delay  t. 
(microseconds) 

2.28 
3 .46 
4.82 
6.25 
7.71 
9.18 
10.6 
12.2 
13.6 


L  =  1 0  cm 

h  =  5  cm 

=  4.422  X  10^  N/m^ 

C.,  =  3.230  X  10^  N/m2 


X 


Angle  of  incidence  incident  SH  wave 

and  angle  of  reflection  of  reflected 

SH  wave  shown  in  y~z  plane  of  (x,y,z) 
coordinate  system. 


Fig.  1 


Transmitting  Receiving 

Transducer  Transducer 


Fig.  3  SH  wave  experiencing  2n-l  reflections  between 

transmitting  transducer  and  receiving  transducer. 


I 

I 
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Fig.  4 


Right  triangle  from  which  trigonometric  functions 
of  0  are  derived. 


9  (degrees) 

Fig.  6  Directivity  function  of  shear  stress  as  function  of 

angle  from  the  z  axis  for  frequency  0.75  MHz. 


function  ^  of  shear  stress  as  function  of 

axis  for  frequency  1.00  MHz- 


as 


Normalized  output  voltage  as  function  of  time  for  f 


.25  MHz. 
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16.  Abstract 

Input-output  characterization  of  fiber  composites  is  studied  theoretically  by 
tracing  SH  waves  in  the  media.  A  fiberglass  epoxy  composite  is  modeled  as  a 
homogeneous  transversely  isotropic  continuum  plate.  The  reflection  of  a  SH  wave 
at  a  stress-free  plane  boundary  in  a  semi-infinite  transversely  isotropic  medium 
is  considered  first.  It  is  found  that  an  incident  SH  wave  reflects  only  a  simi¬ 
lar  SH  wave  back  into  the  medium.  It  is  also  established  that  the  angle  of 
reflection  of  the  reflected  wave  is  equal  to  the  angle  of  incidence  of  the  inci¬ 
dent  wave.  The  phase  velocity  of  the  SH  waves  and  the  delay  time  of  the  SH 
waves  in  reaching  the  receiving  transducer  are  computed  as  functions  of  a  reflec¬ 
tion  index,  defined  as  the  number  of  reflections  of  the  SH  waves  from  the  bottom 
face  of  the  continuum  plate.  The  directivity  function  corresponding  to  the  shear 
stress  associated  with  the  SH  waves  in  the  continuum  plate  is  also  derived  as  a 
function  of  the_ ref lection  index.  A  theoretical  output  voltage  from  the  receiv¬ 
ing  transducer  is  calculated  for  a  tone  burst  (a  periodic  input  voltage  of  finite 
duration).  The  output  voltage  is  shown  for  tone  bursts  of  duration  60  ps  and 
center  frequencies  0.75,  1.00,  and  1.25  MHz.  The  study  enhances  the  quantitative 
and  qualitative  understanding  of  the  nondestructive  evaluation  (NDE)  of  fiber 
composites  which  can  be  modeled  as  transversely  isotropic  media. 
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